Hysteresis and Noise from Electronic Nematicity 
in High Temperature Superconductors 
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An electron nematic is a translationally invariant state which spontaneously breaks the discrete 
rotational symmetry of a host crystal. In a clean square lattice, the electron nematic has two 
preferred orientations, while dopant disorder favors one or the other orientations locally. In this 
way, the electron nematic in a host crystal maps to the random field Ising model (RFIM). Since 
the electron nematic has anisotropic conductivity, we associate each Ising configuration with a 
resistor network, and use what is known about the RFIM to predict new ways to test for electron 
nematicity using noise and hysteresis. In particular, we have uncovered a remarkably robust linear 
relation between the orientational order and the resistance anisotropy which holds over a wide range 
of circumstances. 

PACS numbers: 72.70.+m, 74.25.Fy, 75.10.Nr 



In the high temperature superconductors, in addition 
to superconductivity, there may exist various other types 
of order which break spatial symmetries of the underly- 
ing crystal, especially in the "pseudogap" regime at low 
doping. P, 0. y. 0]. 0| However, it is often surprisingly 
difficult to obtain direct experimental evidence which 
permits one to clearly delineate in which materials, and 
in what range of temperature and doping, these phases 
occur. One such candidate order is the electronic ne- 
matic, which breaks orientational, but not translational, 
symmetry. |3 Orientational long range order (LRO) in- 
duces transport anisotropy, since it is easier to conduct 
along one direction of the electronic nematic than the 
other, and this is a natural way to look for nematic order. 
On the other hand, quenched disorder couples linearly to 
the order parameter, and so, for quasi-2D systems, will 
typically change a true thermodynamic transition into a 
mere crossover. 
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(a) Nematic patches 



(b) Resistor Network 



FIG. 1: Mapping of nematic patches to random resistor net- 
works, (a) Configurations of nematic patches generated by 
the random field Ising model, (b) Corresponding resistor net- 
work, modeling local anisotropic conduction in each nematic 
patch. Solid lines are small resistors, and dotted lines are 
large ones. 
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In this letter, we propose methods for detecting local 
electronic nematic order (nematicity) even in the pres- 
ence of quenched disorder. We show that an electron 
nematic in a disordered square lattice maps to a random 
field Ising model (RFIM) jfj and that the transport prop- 
erties of the system can be determined from a random re- 
sistor network whose local resistances reflect the local ne- 
matic transport anisotropy. In a host crystal such as the 
cuprates, an electron nematic (which may arise from local 
"stripe" correlations) tends to lock to favorable lattice di- 
rections, often either "vertically" or "horizontally" along 
Cu-0 bond directions. The two possible orientations can 
be represented by an Ising pseudospin, a = ±l,[g and 
the tendency of neighboring nematic patches to align cor- 
responds to a ferromagnetic interaction. In any given re- 
gion, disorder due to dopant atoms between the cuprate 
planes produces electric field gradients which locally fa- 



vor one orientation or the other, and act like a random 



field on the electronic nematic, as illustrated in Fig. 1(a) 



Thus an electron nematic in a host crystal maps to the 
RFIM: 
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where J > is the coupling between neighboring ne- 
matic patches. The local disorder field hi is taken to be 
gaussian, with a disorder strength A characterized by the 
width of the gaussian distribution. The symmetry break- 
ing field h may be produced by, e.g., uniaxial strain, high 
current, Q or magnetic field. For example, nematics can 
be aligned by an external magnetic field due to diamag- 
netic anisotropy, in which case h oc o(if 2 ).[j| 

The macroscopic resistance anisotropy of a nematicjlCt 
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(a) R a vs. h. 




(b) -R a vs. ra 



FIG. 2: Hysteresis and comparison of the macroscopic resis- 
tance anisotropy to the orientational order for size L x L = 
100 x 100, disorder strength A = 3J, temperature T = 0, mi- 
croscopic resistance anisotropy r — 2, and external field sweep 
rate f2 = J/N, where N — L x L. (a) Hysteresis of the re- 
sistance anisotropy R a (see text) vs. the symmetry breaking 
field h. (b) Resistance anisotropy R a vs. orientational order 
parameter m. The resistance anisotropy is a remarkably good 
indicator of orientational order. 



Ull ll^ l transforms under rotations in the same way 
as the orientational order parameter to = (^/N)J2i a i 
and is a natural candidate for measuring nematic or- 
der. The (normalized) resistance anisotropy is R a = 

(£t) {.TT-lh)- where r = R *T/Kt is the ratio 
of the extremal macroscopic resistances in the two fully 
oriented states. To obtain the transport properties, we 
map each pattern of local nematic orientations generated 
by Monte Carlo simulations of the RFIM (see Fig. 1(a) I 
to a resistor network (see Fig. 1(b) I which models the 
local anisotropic transport in each nematic patch. Each 
nematic patch (Ising pseudospin) becomes one node in 
the resistor network, with four surrounding resistors de- 
termined by the nematic orientation. For a "vertical" 
nematic patch, we assign the resistors to the "north" 
and "south" to be small, i? sma n, while the resistors to the 
"east" and "west" are large, i?i argo = rR sraa \\. For a "hor- 
izontal" nematic patch, these assignments are reversed. 
When all patches are fully oriented so that m — * 1, 
the macroscopic resistance anisotropy must also satu- 
rate, R a — > 1. More generally, in the thermodynamic 
limit < R a >= mF(m,T,r), where F is an even func- 
tion of to. Remarkably, as we will see below, F = 1 to 
a very good approximation and, under a wide range of 
circumstances, < R a >~ to throughout the entire range 
of m. 

We use a Glauber update method to generate configu- 
rations of the RFIM, with periodic boundary conditions 
on the pseudospin lattice. The details of the algorithm 
are contained in R,ef.[l3j. We then calculate the resis- 
tance anisotropy R a of the corresponding resistor net- 
work by the following method: For R xx , we assign a uni- 
form applied voltage to every site at the far left end of 



the lattice, and a uniform ground to the far right end 
of the lattice, with open voltage boundary conditions in 
the y direction, corresponding to infinite resistors com- 
ing out of the top and bottom edges of the network. We 
then apply the bond propagation algorithm 01 to reduce 
the network to a single resistor equivalent to the macro- 
scopic R xx . R yy is similarly calculated from the same 
starting resistor network, but with boundary conditions 
appropriate for R yy . 

Below a critical disorder strength, the 3D RFIM pos- 
sesses a finite temperature phase transition to a low tem- 
perature ordered phase. However, in 2D the critical 
disorder is zero, and LRO is forbidden. 0] The quasi- 
2D case of coupled RFIM planes also has finite critical 
disorder, ^ but it is exponentially small. We focus here 
on two dimensions, and demonstrate that even in this 
case, where orientational LRO is forbidden, it is possi- 
ble to detect the proximity to order through noise and 
hysteresis measurements. 

We first present results at zero temperature, which ex- 
hibits the nonequilibrium behavior associated with hys- 
teresis. Fig. shows a simulation of the RFIM for a 
system of size L = 100 x 100 at zero temperature. We 
present results for which the field h is incremented at a 
sweep rate = J/N, where N = L x L is the system size 
and h = fit. In Fig. 2(a) we show a hysteresis loop for 
the resistance anisotropy R a vs. the symmetry breaking 
field h, starting from the fully oriented state at h = — oo. 
Fig. 2(b) plots R a vs. to through one cycle of the hys- 
teresis loop. The relationship is remarkably linear, and 
the macroscopic resistance anisotropy R a shows precisely 
the same hysteretic behavior as to. This linear relation 
throughout the cycle is surprising, since for a given mag- 
netization to not too close to 1, there are several values 
possible for the macroscopic resistance anisotropy. This 
is because the flow of current through the random resistor 
network depends not only on the relative concentration of 
"up" and "down" resistor nodes, but also on their spatial 
relation. In fact, for to = 0, R a can take values between 
±(r 2 — l)/(r 2 + 6r + 1). However, for typical configura- 
tions generated by the RFIM, the main contribution to 
R a is controlled by the number of up and down nodes, 
with their spatial relation being a small effect which goes 
to zero as the system size L — > oo or as r — > 1. (We do 
observe variations in R a between configurations with the 
same magnitude of to at finite temperature or for small 
system sizes.) For the same reason, in confined geome- 
tries there can be more noise in the resistance anisotropy 
than is present in the true orientational order parameter. 
The fact that < R a >ss to over a wide range of param- 
eters means that noise and hysteresis in this measurable 
property can be used to detect the electronic nematic. 

Hysteresis subloops can be used to determine qualita- 
tively the relative importance of interactions and disor- 
der. Fig. |31 shows the behavior of subloops in the hys- 
teresis curve of orientational order to vs. orienting field 
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(a) Hysteresis 



(b) Subloops 



FIG. 3: (a) Field ramp-up from a disordered configuration 
at T — 0, with hysteresis subloops, for system size L x L = 
100 x 100 and disorder strength A = 3 J. The arrows indicate 
the order in which the external field sweep is taken. Subloops 
A and B are both taken between the same extremal orienting 
fields h = 1.5J and h — 2.1J, but with different histories, 
(b) Incongruency of subloops A and B, due to the presence of 
interactions. Here the orientational order parameter has been 
shifted to compare the shape of the two subloops. Subloop 
A, executed at lower m, has a higher slope than subloop B, 
executed at higher m. 



h. Starting from zero field and a thermally disordered 
configuration, the field is swept up along path 1, then 
switched back along path 2 to take subloop A before 
continuing along path 3. Then path 4 begins a second 
subloop, within which path 5 begins subloop B, before 
continuing to raise the field through path 6 until m has 
saturated, m — ► 1. 

Notice that the subloops close, and also that once a 
subloop has closed, continuing to raise the field does not 
disturb the structure of the outer loop. This is indicative 
of return-point memory, a characteristic of the RFIM.[l7j 
Additionally, subloops A and B in Fig. |3| are a com- 
parison between the same two field strengths. The fact 

in- 



that the two are incongruent as shown in Fig. 3(b) 



dicates that interactions are present, and possibly even 
avalanches, i.e. that the hysteresis is not simply a linear 
superposition of elementary hysteresis loops of indepen- 
dent grains as in the Preisach model. 0, ll^. Eoj Here, 
we have used a disorder strength of A = 3J > A^ t = 
2.16J. At lower disorder, the subloops become narrower 
and difficult to resolve. Also, as temperature is increased, 
thermal disorder means that subloops become narrower, 
and no longer close precisely. 

Recent magnetization measurements by Panagopoulos 
et al. 0, H2| on L^-zSr^CuC^ reveal hysteretic be- 
havior reminiscent of the RFIM. Although some of the 
results look very similar to our Fig. |3 in fact our model 
does not explicitly include the ferromagnetic moments 
measured in this experiment. However, the small mag- 



nitude of the measured magnetization may be consistent 
with ferromagnetic moments arising at defects in the lo- 
cal striped antiferromagnetic order of a nematic patch. 

We now discuss equilibrium fluctuations. In confined 
geometries such as nanowires and dots, this model can ex- 
hibit telegraph noise in R a due to thermal fluctuations of 
large correlated clusters. Recent transport experiments 
on underdoped YBCO nanowires by Bonetti et al.^j^ re- 
veal telegraph-like noise in the pseudogap regime. Using 
a YBCO nanowire of size 250nmx500nm, they found that 
a time trace of the resistance at constant temperature 
T=100K shows telegraph-like fluctuations of magnitude 
0.25%, on timescales on the order of 50 seconds. These 
large scale, slow switches can be understood within our 
model as thermal fluctuations of a single correlated clus- 
ter of nematic patches. Individual nematic patches ther- 
mally fluctuating between two orientations also produce 
noise, but at higher frequency, and with smaller effect on 
the macroscopic resistance. By comparing the magnitude 
of the smallest typical resistance change to the resistance 
change at a telegraph noise event, one can estimate that 
the correlated cluster in the nanowire contains at least 
5-6 nematic patches. 

The nanowire represents a rather small system size 
when mapped to the RFIM. Neutron scattering exper- 
iments on the incommensurate peaks of underdoped 
YBCO indicate a coherence length of roughly 40nm. If 
we take this as an estimate of the size of one nematic 
patch (mapped to a single pseudospin in the RFIM) , the 
nanowire is about 6 x 12 patches wide, and the effect 
of flipping a correlated cluster of pseudospins can dom- 
inate the response. In Fig. |4(a)[ we show a time series 
of R xx in thermal equilibrium [h = 0) for a small sys- 
tem of size L x L = 6 x 6 at finite temperature with 
disorder strength A = 2 J, which is less than the 3D crit- 
ical value. Notice the sizable thermal fluctuations in the 
macroscopic resistance. The high frequency noise is due 
to the thermal fluctuations of a single Ising pseudospin. 
The lower frequency telegraph noise, in which the resis- 
tance changes dramatically, is due to the correlated fluc- 
tuations of a cluster of pseudospins, in our case a cluster 



of 4 pseudospins. The histogram in Fig. 4(b) shows that 
the system switches between two main states. 

In general, correlated clusters will thermally switch 
with a timescale which rapidly increases with their size. 
Confined geometries such as the nanowire are more likely 
to have a system-spanning cluster which causes slow, 
large scale noise. For larger system sizes, the histogram 
becomes smoother and can have multiple peaks. For sys- 
tems larger than the Imry-Ma correlation length, large 
scale switches are highly unlikely. A probability distri- 
bution of size of cluster vs. frequency can be obtained 
from a histogram of the differential resistance, in which 
large scale switches will always be relegated to the tails 
of the distribution. 

A further way to test for nematic patches is to measure 



4 




(a) R xx vs. time (b) Histogram 



FIG. 4: (a) Telegraph noise in a confined geometry. Resis- 
tance R xx vs. time of a small system of size L x L — 6 x 6 
with disorder strength A = 2J, at temperature T — 0.5J, in 
zero applied field, H — 0. An initial random state is allowed 
to thermalize for 10,000 MC steps before measurements are 
taken. The local microscopic conductivity anisotropy in the 
resistor network is r = 2. Small fluctuations are due to a sin- 
gle nematic patch fluctuating. Large switches are due to the 
thermal fluctuations of a single correlated cluster of nematic 
patches, (b) Histogram of the time series. TV is the number 
of occurrences of a particular resistance. There are two main 
states, one for which the correlated cluster is "up", and the 
other for which it is "down". 



cross correlations between the thermal noise in R xx and 
R yy as a function of time. |24| This can distinguish be- 
tween macroscopic resistance fluctuations due to, e.g., 
fluctuation superconductivity (which would cause cor- 
related fluctuations of R xx and Ryy), or due to RFIM- 
nematic physics (which would cause anticorrelated fluc- 
tuations). 

This general method of connecting local anisotropic 
properties of an electronic nematic to macroscopic behav- 
ior can be extended to many experimental measurements. 
For example, applying an in-plane magnetic field should 
drive 4-fold symmetric incommensurate spin peaks in 
neutron scattering into a 2-fold symmetric pattern, in 
a hysteretic manner as the field is rotated from one Cu- 
O direction to the other. Superfmid density anisotropy 
should also display similar hysteresis with in-plane field 
orientation. Furthermore, the presence of correlated clus- 
ters in the RFIM has implications for STM. Whereas 
small correlated clusters have faster switching dynamics, 
large clusters are much slower. |l8| so that different size 
clusters have different local power spectra. STM can be 
used to do scanning noise spectroscopy, |25l| by measuring 
the power spectrum as a function of position, in order to 
produce a spatial map of correlated nematic clusters. 

In conclusion, we have mapped the electron nematic 
in a host crystal to the random field Ising model. Us- 
ing a further mapping to a random resistor network, we 
have predicted new ways to detect the electron nematic 
in disordered systems. We have demonstrated that the 



macroscopic resistance anisotropy is a good measure of 
orientational order and is expected to display hysteresis 
and thermal noise characteristic of the random field Ising 
model. Recent experiments on noise in YBCO nanowires 
and hysteresis in LSCO exhibit behavior reminiscent of 
this model. 
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